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We calculate a force due to zero-temperature quantum fluctuations on a stationary object in a moving super- 
fluid flow. We model the object by a localized potential varying only in the flow direction and model the flow by 
a three-dimensional weakly interacting Bose-Einstein condensate at zero temperature. We show that this force 
exists for any arbitrarily small flow velocity and discuss the implications for the stability of superfluid flow. 
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Although there are various definitions of superfluidity OL 
one of the defining features of a superfluid is the existence 
of a critical velocity below which the superfluid flows with- 
out dissipation. Landau argued that, by performing a Galilean 
transformation on the ground state of a uniform superfluid, the 
superfluid would become unstable above a well-defined crit- 
ical velocity due to the creation of quasiparticles \2]. If one 
assumes that it is only through the creation of quasiparticles 
that dissipation can occur at T = 0, then one can infer that 
a stationary object in a slow-moving superfluid (with a flow 
velocity well below the critical velocity), would remain in a 
metastable stationary state as there would be no force acting 
on this object. 

In this letter, we show that the phenomenological picture 
of superfluid flow - namely the existence of a metastable state 
below a critical velocity - is incomplete and problematic. We 
illustrate this by examining the case of a localized potential 
fixed in the flow of a three dimensional dilute Bose-Einstein 
condensate. Specifically, we show that a force arises from 
the scattering of zero-temperature quantum fluctuations, an 
Landau ignored in his argument for a critical velocity. We 
demonstrate the existence of this force in an infinitely ex- 
tended condensate at all nonzero flow velocities (including 
velocities much lower than Landau's critical). 

Casimir j3| first showed that zero-temperature quantum 
fluctuations in an electromagnetic (EM) vacuum give rise to 
an attractive force between two closely spaced perfectly con- 
ducting plates. A Casimir-like force, Fbec, can be shown 
to arise from the zero-point quantum fluctuations in a dilute 
BEC, where infinitely thin and infinitely repulsive plates im- 
mersed in a zero-temperature three-dimensional dilute BEC 
replace Casimir's perfect conducting plates. Fbec is given 
by (to leading order) 

FBEC *-iwrf (1) 

where c s is the speed of sound in the dilute BEC, d is the dis- 
tance between the plates, and a is the area of the plates. Note 
that the Casimir and Casimir-like forces in the EM vacuum 
and the BEC vacuum, respectively, have the same expression 
except that the speed of sound replaces the speed of light in 
the BEC case ||4j]. (The expressions also differ by a factor of 
2 due to the fact that EM modes have two transverse polariza- 
tions whereas phonons have only one polarization.) Both of 



these forces arise because boundary conditions are imposed 
on the quantum fluctuations. 

Similarly, we posit that a Casimir-drag force exists on an 
object in a moving dilute BEC. No direct EM analogy can be 
drawn for this situation because no absolute rest frame (where 
the relative motion of an object can be measured) exists for 
an EM vacuum |5[. Nevertheless we maintain that in a su- 
perfluid flow at zero temperature, modeled as a weakly in- 
teracting BEC, around a stationary object (or, by a Galilean 
transformation, a moving object in a stationary superfluid), a 
Casimir-like force should arise due to the boundary conditions 
on the quantum fluctuations of a BEC vacuum. In the specific 
case of a weak potential varying only in the flow direction, we 
show that this drag force exists at all nonzero flow velocities, 
i.e. the effective critical velocity, defined as the velocity below 
which the flow is dissipationless, is zero for this system. 

We now calculate the drag force arising from these quan- 
tum fluctuations. Momentum is not, in general, conserved in 
our system because the stationary object, which is described 
by the potential $(r), breaks the translational symmetry |8]. 
In general, a force on a moving object described by a poten- 
tial $(r) can be written in second quantized notation at zero 
temperature as 



F 



d 3 r(ft (r)[V$(r)]V>(r)) 



T=0: 
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where ip(r) and ^'{r) are field operators that describe the 
weakly interacting BEC flow and obey the standard boson 
commutation relations and the expectation value is taken at 
T = 0. T = is not well defined in the scattering prob- 
lem discussed in this letter so one can view this simply as a 
convenient label of the quantum state that we define in detail 
below. 

We model the superfluid as a weakly interacting three- 
dimensional condensate characterized by an interparticle con- 
tact pseudopotential, gS^^r), where g is determined by the 
2-particle positive scattering length a sc and the mass m of the 
atoms such that g = 4irTi 2 a sc /m. We assume the condensate 
to be dilute such that J pqo? sc <C 1 where po is the condensate 
number density. 

To calculate the force due to quantum fluctuations on a sta- 
tionary object in a superfluid flow, we assume for simplic- 
ity that the object is described by a weak symmetric poten- 
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tial that varies only in the flow direction (which we take as 
the ir-direction), i.e. <£>(r) = r)$>(x) where $(x) = <&(— x) 
and 77 <gc 1 (The parameter 77 gives the order of magnitude of 
the external potential.) This situation, which can, in princi- 
ple, be realized in current dilute BEC experiments, is a spe- 
cific case chosen to show the existence of a finite Casimir-like 
drag force at any arbitrarily small velocity. We will place fur- 
ther restrictions on this potential in the course of this letter as 
needed. 

Because we consider a potential varying only in the flow 
direction, the integrand in eq. (0 is only a function of the 
positional variable x (the y and z dependencies implicit in 
tp(r) and ip ( r ) cancel out), which allows the simplification 
of the force expression to 
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T=0, 
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where A is the cross-sectional area of the object in the flow. 

Although the lack of translational symmetry (due to the 
presence of the object) makes the existence of a drag force 
possible, it does not imply that there will necessarily be a 
drag force. For example, if the small quantum fluctuations 
are ignored, the bosonic field operator tp can be approximated 
by the classical mean field *S?W [9], whose behavior is deter- 
mined by the Gross-Pitaevskii equation (GPE). Working, as 
we will do throughout this letter unless specified otherwise, in 
dimensionless variables in which the length scale is normal- 
ized by the healing length given by (87rpoa sc ) _1 ^ 2 and W is 
normalized by -^/po, the GPE can be written as 

(f + $(sb) -n)V<®(r) + \¥°\r)\ 2 ¥ \r) = 0, (4) 

where T = — V 2 + \/2iqj^ + q 2 /2, the dimensionless speed 
is given by q — c/c s , c is the speed of the flow at x = 00, 
c s = yj pog/m is the speed of sound, and fi = 1 + q 2 is 
the chemical potential (determined by imposing ^°>(r) = 1 
at x = 00). The mean field force arising from the potential 
$(x), given by -Ar\ J dx\^(r)\ 2 ^^ , can be shown to be 
zero below a certain critical flow velocity. This critical flow 
velocity (as measured far from the potential) in a nonuniform 
medium is always smaller than Landau's critical velocity in a 
uniform medium (which in the dilute gas is given by the speed 
of sound) due to nonlinear effects such as vortex shedding [10] 
or the creation of gray solitons 111 ill , and is the velocity at 
which the maximum local fluid velocity reaches the speed of 
sound 1 12]. 

If we go beyond the mean field approximation and take into 
account quantum fluctuations, the bosonic quantum field op- 
erator %jj can be split into a macroscopic classical field "I^ 1 ) 
and a small quantum fluctuation operator <j>: tp = \pW + (p. 

is an improved approximation of the condensate wave 
function as compared to because it includes the effects 
of the quantum fluctuations. Including the fluctuation oper- 
ator in the analysis leads to a depletion of the ground state 
JT3I1 and a correction to the ground state energy irHll . both on 
the order of the diluteness parameter */ poa^ c , which must be 



small in order for the Bogoliubov theory in this paper to be 
valid. The effects of this quantum depletion and their correla- 
tions are well known (see references in |@]). 

The force due to the boundary conditions imposed by the 
potential on the quantum fluctuations is, in general, not zero 
(even below the critical flow velocity in a nonuniform system 
given by the GPE) and can be written as 



F X = -A V J dx(\*W(r)\ 2 +{ft(r)$(r))T=o] 



dx 



(5) 



The fluctuation operator can be expanded in terms of a k 
and al - the quasiparticle annihilation and creation operators, 

respectively - such that cp(r) — J2k — vl(r)ctlj , 

where the sum excludes the condensate mode. For our sys- 
tem of weakly interacting particles to be described by the non- 
interacting quasiparticles, the quasiparticle amplitudes, u k (r) 
and Vk (r), must obey the Bogoliubov-de Gennes (BdG) equa- 
tions lii tn, 



Cu k (r) - {¥V) 2 v k (r) = E k u k (r) 



C*v k {r) - (*W*) 2 u fc (r) - -E k v k (r), 



(6) 



(7) 



where t = f + - /i + 2|* (1) | 2 and k is the dimension- 
less momentum normalized by the healing length. The energy 
eigenvalue for the moving BEC flow is E k = \[2qk x + Eb 
where the Bogoliubov dimensionless dispersion relation for a 
BEC at rest is given by Eb = ky/k 2 + 2. 

Since the BdG equations describe an effective scattering 
problem for the quasiparticle amplitudes, we can solve for 
Mfc(r) and Vk(r) by specifying the incoming u k {r) and v k (r). 
We impose a condition (see II 1 611 for details) on the incoming 
Uk (r) and v k (r) in terms of measurable quantities far from the 
potential that fully determine the quantum state in this prob- 
lem, i.e. the quantum state used in the expectation value in eq. 
0. Despite the fact that the state at zero temperature is tech- 
nically not well defined for this scattering problem, this label 
remains a convenient way to denote the quantum state relevant 
to eq. because the state is annihilated by a. The expecta- 
tion values taken with respect to this state (or at "T = 0") 
can be written in terms of the quasiparticle amplitudes, i.e. 

#(r)0(r)) T=o = E fe kW| 2 - 

The condensate wave function modified by the quantum 
fluctuations is given by the generalized GPE (GGPE) fl7il 

(f + #(«) - /i)* (1) (r) + |* (1) (r)| 2 * (1) (r) - x(V)* (1) (r) 



[2\v k (r)\ 2 *W(r) -u fe (r)^(r)* (1) (r) 



= 0. 



(8) 



where the term proportional to Uk{r)v k (r) is ultraviolet di- 
vergent because of the contact potential approximation and 
must be renormalized (see references in 11810 . The term 
x{r)^^(r) ensures the orthogonality between the excited 
modes and the condensate II 811 and is given by x( r ) — 



3 



EfeCfe^(r) where c k = J d 3 r\^ (r)| 2 OJ«*(rH(r) + 
(r)vl(r)). Because of the properties of $(a;) (described 
below), neither x(r)\I/W(r) nor the renormalization term 
contribute to the dominant order of the calculation below. Be- 
cause we assume the condensate to be dilute, the fluctuation 
terms are small (on the order of J poa 3 c ) so the GGPE (eq. 
[8} can be approximated by the GPE with an effective com- 
plex potential given by ((x) = J2k 2\vk(r)\ 2 - Uk{r)v* k (r) 
HI 911 - (Note the y and z dependencies manifest themselves 
only in the phases of the quantum amplitudes and thus cancel 
out when the phases cancel out). 

To calculate the force arising from zero-point quantum fluc- 
tuations F x , we must first solve the Bogoliubov equations for 
the quantum amplitudes. To obtain a finite force we find it 
convenient to assume /_ dx&(x) — $(0) = where <J>(A) 
is the potential in Fourier space. Extracting the trivial phase 
factor, i.e. Uk(r) = U{x)e b ^ i ' T and Vk{r) — V{x)e l ^ L ' r , and 
working in Fourier space defined by U(x) = J dXe lXx U(X) 
and V(x) = j dXe lXx V(X), the Bogoliubov equations can 
be solved perturbatively to give the quasiparticle amplitudes 
to first order in ry as Ui(k, A) = $(A) sign(^) and 

Vi(k,X) = $(A)^g$sign(i£). T v (k,X) and T v (k,X) 
are quantities easily derived from the Bogoliubov equations 
but, for the sake of clarity, we have chosen not to write out 
their full expressions as they would contribute little to the 
discussion. The sign of the group velocity of the reflected 
quantum fluctuation, denoted by sign(w^), arises from the 
boundary conditions where we exclude exponentially growing 
scattered waves and exclude incoming scattered waves due to 
causality. The group velocity of the reflected wave is given 
by v R = Isfk where the wavenumber of the reflected mode 

J 9 oar 

is given by Xr = k x + kji and is given by the non-trivial 
real root of the characteristic equation of the coupled Bogoli- 
ubov equations C(k, A) = A [A 3 + Ak x X 2 + (4fc 2 + 2k 2 + 2- 
2q 2 )X + Ak x (k 2 + 1) + 2^2qE B \ = 0. Assuming k x = kf, 
then vf > if -1 < / < f c and vf < if f c < f < 1 

where f c = - %£g±L 

Next, these quasiparticle amplitudes can be used to deter- 
mine the effective complex potential ((x) in the GGPE to give 
&W(r). Then, integrating over all momenta of the quantum 
fluctuations, the Casimir-like force due to the quantum fluctu- 
ations F x can be divided into two contributions as seen in eq. 
(0, given at the dominant order in rj as 



F x = - d 3 k(F cond (k) + F fluc (k)) 



(9) 



where F x — F x /rj 2 poAyJ poa 3 c and the zeroth order interac- 
tion pressure is given by po = gp^/2. The contribution to the 
force due to the condensate modified by the quantum fluctua- 
tions is given by 

F cond (k) = Res x=kn ^$(A)A/(A 2 + 2 - 2q 2 ) 
{[U (k)(l + ^) - 4Vo(fc)]^i(A, fc) + V (k)(l - ^ 
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FIG. 1: The scaled force due to quantum fluctuations, F x = 
F x /r] 2 poA\/poa'i c , acting on a stationary localized potential de- 
scribed by the parameters fco and A as a function of the mach number 
of the flow, q (measured far from the potential, i.e. at x = oo), for 
three different values of ko: fco = 5.0 is given by the solid line, 
fco = 2.4 is given by the dashed-dotted line, fco = 1.8 is given by the 
dashed line. The width of the potential in real (and Fourier) space is 
equal to the healing length, i.e. A = 1. (fco and A are both normal- 
ized by the healing length.) 



and the contribution to the force given directly by the quantum 
fluctuations is given by 



Ffi uc (k) = Res 



^$(A)Vo(fc)AT/!(A,fc), (11) 



where Res^/^ is the residue at A = &r. The zeroth order 
quantum amplitudes (for a homogeneous gas at rest) are given 



by U (k) 



k 2 + l 



1) andV (k) 



fc 2 +i 



1 



E B 1 *-) '"V'V — \ l 2 \ E t 

Finally, to illustrate the calculation of F x , we shall define a 
specific potential describing the stationary object in the flow 
as A(A) = [$(A)] 2 = ifor|A + fc | < A/2, |A-fe | < A/2 
and zero otherwise. In real space this can be written as 
A(x) = sin ^A /2) 2 cos ( fc o^) where A is positive, 1/A is a 
measure of the width of the potential in real space and fco is 
a measure of the typical wavenumber of the potential in real 
space. Both A and fco are normalized by the healing length. 
We assume fco > A/2 so that the potential satisfies the con- 
dition $(0) = 0. Note that as A — ► the potential becomes 
periodic in real space and delocalized; our analysis would then 
no longer apply. F x peaks when the width in Fourier space (or 
real space) is on the order of the healing length, i.e. Awl 
and disappears in the localized and delocalized limits. 

In the rest frame of $(r), the grand canonical energy de- 
creases with increasing flow speed. It follows from this that 
the flow would accelerate in the presence of dissipation, im- 
plying a negative drag force similar to the behavior of a mov- 
ing gray soliton [ 20] . At larger fco such behavior is exhibited 
by F x as seen in Figure ^ At smaller fco or, equivalently, at 
larger characteristic wavelengths of the potential in real space, 
an instability arises as the flow speed increases. This insta- 
bility is consistent with the dynamical modulation instability 
_ observed for a perfect lattice in the GPE 1 21 ], which also oc- 
(Asrtw long characteristic wavelengths (small fco) and large 
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speeds. It is also perhaps instructive to recall the non-trivial 
and highly geometry-dependent sign of the Casimir force in 
the EM vacuum (for example, a set-up of parallel conductors 
in an EM vacuum leads to an attractive force while a cubical 
cavity leads to a repulsive force). 

Even though this Casimir-like force exists for speeds much 
lower than the speed of sound, which in Bogoliubov's theory 
is equal to Landau's critical velocity for the onset of dissipa- 
tion, this drag force does not explicitly violate the spirit of 
Landau's principle Let us recall that Landau's princi- 

ple states that above a critical velocity a system can lower 
its energy by the creation of quasiparticles [2]. In our anal- 
ysis, however, we do not assume that any quasiparticles are 
created, i.e. the system remains at T = in the sense de- 
fined above. This drag force arises from the scattering of 
zero-temperature quantum fluctuations and is not caused by 
the creation of quasiparticles, which must satisfy the Landau 
criterion; instead it is caused by the changing nature of the 
eigenstates of the quantum fluctuations, akin to the original 
Casimir force between two conductors. 

In this letter, we have shown in a specific example that a fi- 
nite Casimir-like drag force arises in a moving BEC at T = 0. 
In this case, unlike for the nucleation of vortices, there is no 
free energy barrier to cross and, at least for this particular sit- 
uation, the effective critical velocity is zero. Since a nonzero 
effective critical velocity does exist at the dominant order (on 
the order of the mean field), one would expect to find the 
semblance of a nonzero critical velocity as seen in 1 23], even 
though, at least in the case considered here, the actual critical 
velocity for the system might be zero. 

We expect this Casimir-like force to act upon any density 
perturbation, including those created by laser fields l23ll . un- 
trapped impurities II24I1 . vortices, etc., moving relative to the 
condensate. We also expect this force to be more apparent in 
condensates of lower dimensions due to the enhancement of 
quantum fluctuations. Finally, although the present analysis 
assumes a dilute gas and does not strictly apply to dense sys- 
tems, we do expect a Casimir-like force from quantum fluc- 
tuations also to exist for superfluid liquid Helium. In fact, 
the Casimir-like force might have a stronger effect in liquid 
Helium since quantum fluctuations dominate the Helium con- 
densate at T = 0. 

We conclude by noting that although we have discussed a 
force that exists on a stationary object in a superfluid moving 
at any arbitrarily small velocity, these results are not incon- 
sistent with the existence of persistent superfluid currents in 
toroidal geometries. In this letter where we consider an in- 
finite medium, the Casimir-like force arises from the nonlo- 
cal perturbation of the scattered quantum fluctuations. These 
scattered fluctuations can be seen to transport energy far from 
the potential similar to wave-drag situations in classical flu- 
ids. However, in a finite geometry such as a superflow in a 
torus the scattered waves will interact with the localized ob- 
ject. In the steady state, we expect these backscattered waves 
to eventually cancel out the effect discussed in this paper and 
thus remain consistent with a persistent superflow. In this sys- 



tem, the experimental manifestation of this Casimir-like effect 
would not be a drag force at arbitrarily low speeds but rather 
the presence of these scattered waves which could, in princi- 
ple, be observed. 
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